Stationary quantum statistics of a non-Markovian atom laser 
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We present a steady state analysis of a quantum-mechanical model of an atom laser. A single-mode atomic 
trap coupled to a continuum of external modes is driven by a saturable pumping mechanism. In the dilute 
flux regime, where atom-atom interactions are negligible in the output, we have been able to solve this model 
without making the Born-Markov approximation. The more exact treatment has a different effective damping 
rate and occupation of the lasing mode, as well as a shifted frequency and linewidth of the output. We examine 
gravitational damping numerically, finding linewidths and frequency shifts for a range of pumping rates. We 
treat mean field damping analytically, finding a memory function for the Thomas-Fermi regime. The occupation 
and linewidth are found to have a nonlinear scaling behavior which has implications for the stability of atom 
lasers. 

PACS numbers: 03.75.Fi,03.75-b,03.75.Be 



INTRODUCTION 

The development of Bose-Einstein condensation (BEC) []lj] 
has opened up the study of the atom laser: a coherent source 
of atoms analogous to the optical laser [0]. These atomic 
sources may demonstrate many of the advantageous proper- 
ties of optical lasers, such as spatial and temporal coherence, 
and high spectral flux. This paper examines a quantum me- 
chanical model of an atom laser consisting of a single-mode 
BEC coupled to a continuous output field, which may be a 
valid model for stable lasing systems when they have reached 
steady state. We find a compact form for the standard quantum 
limit for such a laser model which can be determined with- 
out making the Born-Markov approximation. This leads to 
nonlinear scaling of the steady-state occupation of the lasing 
mode with pumping rate, which has implications for the de- 
sign of linewidth narrowing quantum feedback schemes and 
for the stability of atom lasers. 

Outcoupling from a BEC has been achieved using radio 
frequency (rf) radiation @ and Raman transitions [||] to 
change the internal state of the atoms to a non-trapped state. 
Coupling the atoms out more slowly reduces the linewidth of 
the output at the expense of reducing the beam flux ^ . Op- 
tical lasers achieve a narrow linewidth through a competition 
between a saturable pumping mechanism and the damping of 
the lasing mode. This allows a higher pumping rate to in- 
crease both the total flux and the spectral density of the out- 
put. An atom laser with gain narrowing must also have a sat- 
urable pumping mechanism that operates at the same time as 
the damping [^]. A recent experiment managed to combine 
a continuously produced series of condensates to maintain a 
BEC in a trap indefinitely [ph. This procedure was an excel- 
lent first step towards providing a continuous pumping mech- 
anism. The final step will be to produce a stimulated transition 
into the final BEC, which will allow mode-locking and reduce 



the phase diffusion of the lasing mode. There have been sev- 
eral proposals for producing this kind of transition [gfl, but 
they have not yet been experimentally achieved. 

The analogy between optical lasers and atom lasers has 
some limits. Theoretical examination of the properties of the 
output of atom lasers is complicated by the slow dispersion 
of atoms and the high interactions in atomic systems, which 
means that in general they exhibit multimode behaviour and 
the Born-Markov approximation cannot be made [|l0 , 11 , j~2|] . 
Atom laser models have made either a single-mode approx- 
imation for the lasing state or the mean-field approximation 
[jl]]. The first class of models cannot examine the multimode 
behaviour of the laser, which has been shown to have signif- 
icant effects on the stability of the device jl2[]. The second 
class cannot examine the linewidth of an atom laser, which is 
a function of the quantum statistics of the lasing mode JDh,, 
Many have also made the Born-Markov approximation, which 
has been shown to be invalid. This paper examines the output 
properties of an atom laser which has reached steady state, 
which means that the quantum statistics of the lasing mode 
must be included in the model. In order to avoid a multimode 
quantum field analysis, we assume that the steady state BEC 
can be modelled by a single mode, although we make no as- 
sumptions in our formalism about the spatial structure of that 
mode. 



THE MODEL 

The atom laser is modelled by separating it into three parts. 
Recent semiclassical modelling has shown that multimode be- 
haviour of the BEC can be very important, often determining 
the stability of the system [12|. Once conditions for stabil- 



ity have been achieved, by definition we expect the system to 
converge to a stationary state. While the details of this sta- 
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tionary state will depend on the trapping conditions, pumping 
method and damping rate, we will assume in this work that in 
the steady state the BEC can be adequately described by a sin- 
gle mode. This single-mode BEC, also called the 'lasing' or 
'system' mode, is modelled by the annihilation(creation) op- 
erator at" and Hamiltonian H s . Since the external field and 
trapped atoms are in a different electronic state, the free atoms 
are not necessarily affected by the trapping potential. We de- 
scribe the continuum of external modes with field operators 
Cp ' and Hamiltonian H a . The label p is a parameterization 
of the output eigenstates which in general may be degenerate. 
The coupling between the lasing mode and the output modes 
will be described by the interaction Hamiltonian Hi. At this 
stage, we will describe the pump which irreversibly couples 
the atoms from the pump reservoir into the trap mode by H p . 
The total Hamiltonian is then written 



where 

H s 
Hi 



H tot = H p + H s + Hi + H Q (1) 
Tius a) a, (2) 

(4) 



(5) 
(6) 



H = J dp huj p cpCp , 

and the operators satisfy 

[a,a f ] = 



[c P ,c],,] = 5{p-p) 



with all other commutators zero. The natural description of 
the output field in the dilute regime is in the basis of energy 
eigenstates of the output Hamiltonian. The output field eigen- 
functions satisfy 



TlUJr 



(7) 



In the position representation the effective coupling function 
between the trapped mode and the output field is the prod- 
uct of the condensate wave function and the spatial profile of 
the coupling field. We denote the effective outcoupled atomic 
state ket by \k, t) = J dp k(j>, t)\u p ). We have taken out the 
overall coupling rate yfiy in the interaction Hamiltonian so that 
n(p,t) satisfies J dp \n(p, t)\ 2 — 1, The system is described 
in Fig. U 

To carry out a stationary analysis we assume the coupling 
amplitude is constant in time, apart from the Rabi frequency 
of the coupling which may be eliminated by moving into a 
rotating frame. This corresponds to a constant system wave- 
function and spatial profile of the coupling. 



Reversible output coupling 

Atomic output coupling has been achieved through Raman 
transitions and radio frequency coupling of the trapped atoms 



lasing mode 




FIG. 1 : Schematic of the atom laser model. 

to an un-trapped or anti-trapped magnetic sublevel §]. 
Later, we will consider the rf-coupling since in practice the 
energy shift due to rf-coupling is negligible, unlike the Raman 
process which imparts a significant impulse during coupling. 
Such an impulse is modelled in our formalism by shifting the 
center of coupling function k in momentum space. Any re- 
alistic model of an atom laser should describe the coupling 
process in full as the usual Born and Markov approximations 
which are so useful in quantum optics are not generally appli- 
cable to atomic output couplers [[To|]. In making the Markov 
approximation for the atomic output coupler, one effectively 
assumes that the output process is irreversible, so that outcou- 
pled atoms have only a minor influence on the system dynam- 
ics. While this works well in the optical case, the inertia of 
atoms alters the coupling process in a crucial way. An atom 
may be coupled into an untrapped internal state, but because 
it will linger in the interaction region for some time, and it 
may evolve for a short time in the output potential, and then 
undergo a transition back into the trapped state. The slow 
dissipation of atoms may cause the occupation of the system 
mode to increase. 

We wish to solve the equations of motion of this model in 
the stationary long time limit. To clarify matters we will first 
reiterate some elementary results. The second-quantized out- 
put field operator for the output field can be expanded over the 
orthogonal set of modes 



dp (x\up)c p . 



(8) 



Transforming into an interaction picture leaving only Hi sim- 
plifies the description of the output field. The interaction pic- 
ture operators are 

M*>t) = e< H '<- t -*>V R j>{x,t )e- iH '<- t - t 'V R , 

(9) 

oj(t) = e i(H,+H p )(t-t a )/n a ^ o)e -i(H e+ H p )(t-t yn_ 

(10) 

We set t Q = in the rest of this work. The Green's function 
propagator for the output field is 

G(x,y;t) = [fa(x,t),$(x,t)] - (x\e- iH ^ h \y), 

(11) 

where 

(12) 



tpi(x,t) = dp (x\u p ) c p (0)e" 
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Equations of motion without pumping 
The memory function for the output beam is defined by [ |l7| ] 



F(x,t) 



\k) = / dyG(x,y;t) K (y). (13) 



This may be interpreted as the output wavefunction, after 
evolving for time t in the output field. The Heisenberg equa- 
tions of motion for the output and trap operators are 



dt 

tfat(f) 
dt 



iuj p c£ (*) ~ VjK(p)*a^(t) 



--[a){t),H p ]+iuj a){t) 
/ dp n(p)cl(t). 



(14) 



(15) 



Using the formal solution of ( |l4| ) with the Heisenberg picture 
version of (|J), leads to the linear integral equation for the out- 
put field 



ip(x, t) = ^(x, t) - V7 / dt' a{t')F(x, t - if). 
Jo 



(16) 



This shows how the output field may be constructed if the 
evolution of the trap operator is known. This simple form for 
the output field arises because we are considering the case of 
a dilute output beam. 

The output coupling is reversible so we require a descrip- 
tion of the way atoms are coupled back into the trap. The over- 
lap of the output wavefunction with the trap mode describes 
the probability of this occuring, and is defined by [O] 



/(*) 



,-iH a t/h 



)k\ (17) 



= dx dyG(x,y;t)K*(x) K (y). (18) 



The Greens function for the output potential of H Q and the 
form of the coupling completely describes the influence of 
the output on the system evolution. In the limit where f(t) 
becomes a delta function, and the resulting equations of mo- 
tion will be effectively memoryless or Markovian. This cor- 
responds to the quantum optical limit of the coupling where 
photon loss from an optical cavity may be well described by a 
very narrow memory function 



'ity n 
&■ 



have managed to merge independently produced condensates 
in a pulsed, but quasi-continuous manner, which provides an 
excellent starting point but does not involve the stimulated 
transition [|8]]. The stimulated nature of the transition is nec- 
essary for two reasons. The phase of the lasing mode must be 
preserved, and the small linewidth of a laser is a product of 
the competition between a damping process and a saturable, 
Bose-enhanced pumping [Q]. It is a highly non-trivial experi- 
mental problem to design and build such a pumping process, 
and rather than address these issues, we will use a model pro- 
cess that should exhibit the qualitatively similar features to 
any successful continuous pumping process. A pumping pro- 
cess which can satisfy these requirements can be found in a 
model where cooled atoms in an excited state are passed over 
the trap containing the lasing mode [B], The photon emission 
of the atoms would be stimulated by the presence of the highly 
occupied ground state, and they will make a transition into that 
state and emit a photon Jl5|]. For a sufficiently optically thin 
sample, which can be made possible by having a very tight, 
effectively low dimensional trap, the photon is unlikely to be 
reabsorbed, and the process is effectively irreversible. 

If the pump reservoir is sufficiently isolated from the cavity 
and external fields, and the pumping process is designed to be 
irreversible, then we may trace over the pump reservoir states 
to produce a master equation term for the reduced density ma- 
trix which describes only the cavity and external fields. 

We choose to model an optical cooling process rather than 
the more experimentally successful evaporative cooling pro- 
cess as we are particularly interested in designing a contin- 
uously pumped system with a steady state. The reader is 
referred to the literature for the details of such proposals 
After we have traced over the pump modes, a 
term due to the effect of the pump in the equations of mo- 
tion for (a'(t)) (a'a)(t) is obtained by tracing over the trap 
modes in the number state basis JIH]. It is important to note 
that we are not tracing over the output field modes, so the re- 
duced density matrix spans the output field as well as the trap 
field, The term due to the pump in the master equation for the 
density matrix is 



(P) 



I'll III p 



= rV[a)](n s + A[a^}) 



(19) 



Continuous pumping 



where r is the saturated gain rate and n s is the saturation bo- 
son number. The superoperators are 

A c ] = ^( ctc P-P ct c), 

T>[c]p = cpc 1 — Ap. 

This is the same pumping term as derived by Wiseman JT^]. 
In the trap number state basis we denote p nm — (n\p\m), and 
the pump term becomes 



The biggest hurdle in the development of a mode-locked, 
gain-narrowed atom laser is the production of a continuous 
pumping mechanism which can repopulate the lasing mode 
by a stimulated transition. The best experiments at this stage 



i.Pn.m)% 



n s + (n + m)/2 



Pn — l,m— 1 



(n + m + 2)/2 
n s + (n + m + 2)/2 



Pn 



(20) 
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This leads to the term in the (aJa) equation 



dt 



00 1 1 

En + 1 
I - ! 1 Pn,n 
n + 1 + 



— ' V n + n s (n + n s )(n + 1 + n s ) 

n— v 



JV 



2V- 



(21) 



The last line is obtained by assuming that when the mean 
trapped mode occupation TV becomes large, we may treat n 
as a continuous variable and expand the brackets in a power 
series about N. In contributing to the linewidth it will be seen 
that the correction is of order r/JV 3 since a further division 
by N occurs in deriving the spectrum. Since we consider the 
regime n s < N and 1 C JV we neglect the correction here- 
after. The equation of motion for (af) can be treated similarly, 
but in addition it is assumed that as the system is driven further 
above threshold gain narrowing will occur, permitting the co- 
herent state approximation p ~ |V r /V)(v / JV| [pp. The matrix 
elements become 

Pn,m - 7== H(m| (22) 

Vn!m! 

and the equation of motion for (a'(t)) has the pump contribu- 
tion 



— 

dt ' 



,t( t )) 



pump 



^ 2(n + n s ) + 1 

71=0 V / 



2(JV- 



(ot(t)) = P(ot(t)). 



(23) 



Expanding the exact expression in a power series about 1/ JV 
and using the coherent state approximation shows the correc- 
tion is of order r/JV 3 . 

Steady-state occupation in the Born approximation 

In order to easily compare our results with the familiar 
quantum optical case it will be necessary to discuss the Born 
approximation. The Born approximation assumes that the 
density matrix p can be written as a tensor product of the re- 
duced density matirx for system a, and the reduced density 
matrix for the the output field. It further assumes that the out- 
put field remains in its original state. If we trace over the 
output field modes, we obtain a master equation for the las- 
ing mode which has the same pumping term as given above. 
The reduced density now describes only the cavity mode, and 
o'run = (n\a\m) are c numbers. The damping term has the 
form 

(^)damp = -7 / du {f(u)e luJot [a)ao-{t - u) 
Jo 

-aa(t- u)a f ]+ H.c.}. (24) 



Combining the pumping and damping terms we can find the 
equation of motion for the steady-state trap occupation distri- 
bution,^ = 0- n ,n, 

. . , rn r(n + 1) , . 

Pn{t) = ; Pn-l(t) ; — 7Pn{t) 

n + n s n + n s + 1 



-2n / dwRc(/(u)e w )p„(t-u) 
Jo 

+2(71 + 1) / duRe(f(u)e luJot )p n+1 (t - u). 
Jo 



(25) 



By setting the derivatives to zero and assuming that the func- 
tions p n (t) approach a constant in the long time limit, we 
find the recursion relation 



rn 



-p s n s -i + (n + 1)2 7b ^+i 



r(n + 1) 
n + 71= + 1 



2n lBM p s n s = 0, (26) 



where 



lBM = 7 / du Re(f(u)e lu " u ), 
Jo 

is the effective damping constant. The solution 

■(r/2 7 B„) n 



(n + n s )\ 



(27) 



(28) 



is almost identical to that obtained for the optical laser Jl3|]. 
The distribution is thermal for rf^BM < n s , and for 
r/2j BM ^> n s , it approaches a Poissonian distribution with 
mean N and variance V given by 



27b m 
V = N + n s 



(29) 
(30) 



Frequency shift and linewidth in the Markov approximation 

If the evolution of the system mode is much slower than the 
timescale of decay for the memory function f(t), the Markov 
approximation may be used to obtain a useful limit of the 
equations of motion. Using the solution for c p (t), we may 
write 



da\t) 



dt 



— 7 / du a' (u) f* (t — u) , 
Jo 



where 



C(t) = / dp K{p)*c p {Q)e-^ 



(31) 



(32) 
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is an operator-valued noise term generated by the initial state and 
of the output which serves to preserve the commutation rela- 
tions. The commutator is given by the memory function found 
above dt 



(t) = (P - "Ibm + i(u + A M )) g (1) (i) 



(38) 



(33) 



which decays on a timescale called the memory time Jl7|], 
which we denote by T m . This characterizes the timescale 
for which an output atom may influence the system evolution. 
The Markov approximation is made by moving to a rotating 
frame at the system frequency a(t) = a(t)e~ luJot , c p (t) = 



in which, if T m is very small compared with the 



timescale for evolution of the system operator, we may take 



at 



(t) out of the integrals and obtain the equation of motion 



dt 



[a*, H p ] +iA M a*(t) 



(34) 



-Tbm^W + V^), (35) 



where 



A, 



7 / dulm(f(u)e luJ ° u ). 



(36) 



We can now construct equations of motion for the occu- 
pation and the first order correlation function g^(r) = 
limt_» o(a'(i + r)a(t))/N. In the Markov approximation 
the equations of motion for the occupation and the two time 
correlation are uncoupled, and the equation of motion for 
{a'(t + r)o(t)) is independent of t. The equations are 



dN(t) 
dt 



N(t) 



N(t) + n s 



-2 lBM N(t), 



(37) 



The steady state two-time correlation is an exponential lead- 
ing to a Lorentzian spectrum. Using the notation g^\t) = 
cxp(— F M \t\ +i(A M + uj )t), the steady state parameters are 



N 



r M = 



27, 



2(N + n s ) 2(7V + n s ) + 1 
+ 0(N- 3 ) , 



4(N + n s y 



(39) 



(40) 



and the frequency shift is given by ([36|). ^From this point of 
view it is clear that the Born and Markov approximations are 
closely related if we are only describing the trap; in particular, 
the effective damping constants are identical. This is because 
the output and trap spaces remain uncorrected in both approx- 
imations. 

EQUATIONS OF MOTION 

We are now in a position to find the non-Markovian equa- 
tions of motion including the effect of the pump. The equa- 
tions of motion for the system+output, ( |l5| ) and (|l4|), are lin- 
ear and Markovian; however the equations of motion for the 
system mode alone will depend on the output field dynam- 
ics through the memory function f(t). Using the pump terms 
from Equations (pb and (|23|), and substituting the formal so- 
lution of ( p"4] ) into (|15|), the equations of motion for the trapped 
mode are 



4(aM(i) = r-J^-_- 7 2Rc( f du f{u)(a\t)a(t - u)) 
dt N + n s \J 

^-(a\t + T )a{t)) = {iuj + P){a\t + T)a{t)) - 1 ( + duf*(t + t - u)(<J (u)a(t)>, 
or J 



(41) 
(42) 



and for the output modes we find 



d(4c p ) 
dt 



7 |k(p)| 2 2Rc( J due-*"*(*- u >(a + (t)a(u)) 



(43) 



These equations of motion were first derived in JTT|], with a 
seemingly minor modification of the pump term in (|4l|). This 
will be discussed in the following section as it turns out to 
have important implications for the predictions of the model. 
The first two coupled integro-differential equations describe 
the system and output dynamics, starting from when the pump 



is initiated. We are primarily interested in the stationary so- 
lutions, and once the system reaches a steady state the output 
energy spectrum will not directly be of interest as it is steadily 
growing. The long time limit of ( p3[ ) is the spectrum of the 
output energy flux which will become stationary. In this case 
we will see that the relationship between the output flux spec- 
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trum and the power spectrum of the trapped mode is in keep- 
ing with Fermi's golden rule. 

The equations of motion (37) and (38) of [flip are identical 



to our equations (|41[ ) and Q42| ) apart from the additional factor 

N/(N- 



■n s ) in (|41|). The treatment of the pump term in [|1J 
where this factor was taken as unity, produced incorrect results 
because a similar approximation was not made in the pump 
term in (38) of [Q. We will discuss this inconsistency in 
more detail once we have derived the stationary equations of 
motion. 



Stationary equations of motion 



The stationary solutions of ( p4| ) and (^2|) may be found by 
breaking up the integrals and using the fact that the memory 
function vanishes in the long time limit. We will assume the 
system reaches a steady state at t = t 3S . The integral in ( filj ) 
may then be written as 



lim 

t — >oc 



du /(it) (aJ (t)a(t — u)) = lim 



du f(u) g {1 \u)N + lim 

t— >CO 



du f(u) g {1 \u)N. 



du f(u) (a^(t)a(t-u)) (44) 



(45) 



o 



Similar manipulations for the integral of (42 ) lead to the term 

rt+T 



lim 



du f*{t + r-u) {a\t)a{t-u)) = I du f(u)* gW (t - u)N, 



(46) 



where we have assumed stationarity when t is much larger than the memory time T m , defined as the support of the memory 
function. Using these expression we find the stationary equations of motion for the trapped mode as 



( J duf(u)gW(u) ) , 



dr< 



,U) 



(r) - (iw + P)g w (T)-7 / du/*(u)gW(T-u), 





dt 



= 2^ 7 iV|( Wp |Ac)| 2 S( Wp ). 



(47) 
(48) 
(49) 



which are exact in the long time limit. Equations (|47|) and 
( |48| ) may be solved self consistently for the mean occupa- 
tion and the two time correlation of the trap mode. The vari- 
ables are the pumping rate r, the saturation Boson number 
n s , the damping rate 7, the memory function f(t), and the 
system frequency u . One might suspect that the equation 
for g^ 1 ' (rg^ ( would be amenable to some kind of transform 
method; however, it is easily seen that this approach does not 
lead to any major simplifications. If the integral is decom- 
posed into (J Q + J T ) dr, the first term is a convolution, but 
the second term is still important for small r, and this is the 
region of interest. The present form has the added simplic- 
ity that we only require the memory function f(t), and not its 
Laplace transform which is generally more complicated. 

In deriving the output flux expression we have defined the 
stationary power spectrum of the trapped mode as 



S(w) 



lim 

t — >oo 



4 (aHt + r)a(t)) _. 
get — - — , 1 , — — e 

-t 



(50) 



^/drgW(r) 



(51) 



Equation (|49|) is then readily found from (J4-3J) by handling the 
integrals in a similar manner as for the trapped mode equa- 
tions. This expression is essentially Fermi's Golden rule for 
the atom laser model. It is easily verified that this result im- 
plies that our stationary results are exact within first order 
perturbation theory [[l8||, This is because the coupling only 
has non-zero matrix elements between the trap and output 
eigenspaces, so that all higher order perturbation expressions, 
which arise from energy nonconserving virtual transitions be- 
tween output eigenstates, vanish identically. 

If the spectrum for the trap is an exponential g^(r) = 
exp (— f |r| + ilot), corresponding to a Lorentzian spectrum, 
then in the limit that the system linewidth becomes very nar- 
row the flux spectrum becomes 



d{c\c p 
dt 



2ttjN\(u p \k)\ 2 S((jJ p ~ u). 



(52) 
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The amplitude of the output spectral flux is now determined 
by the matrix elements | (u p \k}\ 2 . Using ( p^ ) we may find the 
corresponding expression for the spatial correlation 

|<vV)v>(*)> 



= 2n 1 N J dp<t>* p (x')<t> p (x)\(up\K)rS(u P ), (53) 

where <fi p (x) = (x\u p ) is the output eigenfunction in the po- 
sition representation. In deriving these results the interactions 
in the output beam and the trap have been neglected; however, 
because the equation for c p (t) remains linear in a(t) when 
the interactions for the trapped mode are included, Equations 
( p^ ) and (|3j) for the output flux also hold when interactions 
are significant for the trapped mode. The beam may be con- 
sidered a true atom laser beam when the spatial coherence 
length given by ( |53| ) is significantly longer than the thermal 
de-Broglie wavelength of the atoms [[7]]. It is clear from ( |53"| ) 
that this will place restrictions on the spectral width of the trap 
and the force that acts on the output atoms. As the spectrum 
approaches a delta function in frequency, the coherence length 
will be determined by the eigenfunctions that have the same 
energy as the trapped mode. 

Since we are interested in a highly occupied trapped mode 
our description of the trap should, in the single mode model, 
include a self-interaction term of the form H co \\ = hCa^a'aa 
in the Hamiltonian, where C is the nonlinear interaction 
strength. In the limit where the interactions dominate, the 
spectrum becomes 



g ( 1 )(r) = exp[-iV(l 



2iCr 



)]• 



(54) 



which has the interesting feature of periodic revivals at r = 
mir/C, for integer m fll9[]. In this regime the non-Markovian 
term in ( |42[ ) will become unimportant. However, there are 
regimes of interest where the interactions are weak, corre- 
sponding to a very weak harmonic confinement, or to rela- 
tively low occupation number. In this regime, the Markovian 
effects will alter the effective occupation, producing signifi- 
cant linewidth and frequency changes in the spectrum. 



The Self-Consistent Markov Approximation 

Although the traditional Born-Markov approximation can- 
not be made for atom lasers, as the linewidth becomes narrow 
there is still a timescale separation between the memory time 
of the outcoupling process and the temporal coherence of the 
lasing mode. This allows us to deal with the effects of the sys- 
tem memory using the self-consistent Markov approximation, 
which was first presented in [^TJ. This involves assuming the 
spectrum is Lorentzian and solving for the linewidth and spec- 
tral shift self consistently in an appropriate rotating frame. We 
are now in position to find the sufficient condition for the va- 
lidity of this approximation in the steady state. Since the form 
of gt 1 ' (t) in our equations of motion is still unknown, we must 



make some assumption about its form, if we are to progress 
further. However, since we are interested in the regime where 
linewidth narrowing is expected to occur, we will assume that 
the memory function f(t) has a short decay time compared 
with the decay of g^'(t). In this regime the decay of g^(t) 



may be neglected in the integrals of (WTb and (48 ). This corre- 



sponds to the case where the atoms are ejected from the trap 
rapidly compared with the system coherence time. Making the 
ansatz g^'(t) = cxp (— r SJV/ |i|) exp (i{w + A SM )t), where 
T SM and A SM are real constants, the equations of motion be- 
come 



r = (N + n s )2-y 



du Rc 



(/(«)< 



i(u +Asm)u 



(55) 



poo 

-7/ du f*{u)e- l(u} " +AsM)u e VsM{t -^- u]) . (56) 



As T~lj becomes smaller than T m we can ignore the vari- 
ation of the rightmost factor in the second integral. This ap- 
proximation explicitly removes any time dependence from the 
equation, as required for the consistence of our ansatz. If we 
define the following pair of self-consistent parameters that de- 
pend solely on the memory function: 



Ism = 7^ du Re (/( u )e^°+ As -)«) (57) 
A SAJ = 7 j du Im (/( u )e^° +AsM > M ) , (58) 



then we can immediately determine the steady-state lasing 
mode occupation and linewidth of the atom laser: 



N 



27s* 



2(N + n s ) 2(N + n s ) + l 
0(N~ 3 ). 



A(N + n s f 



(59) 



(60) 



These equations now take a very similar form to the Markov 
equations, but must be solved self consistently for the system 
variables in the steady state. The possibility of a frequency 
shift also arises when the more severe Markov approximation 
is made, but this only depends on the bare trap frequency, and 
corresponds to the limit A SAf <C lo of (|8|). (Note that the 
linewidth, while taking the same functional form as the usual 
Born-Markov result, depends on N which is determined by 
(|5^), and hence will be changed when A SM 7^ A M ). 

Equations ( |58| ) and ( |60| ) are our main results. Comparison 
with the Markov expressions shows that the frequency shift 
for the system mode and the effective damping constant must 
now be found self consistently. If the output damping is dom- 
inated by gravity and the wavefunction of the system state is 
weakly dependent on N the memory function will not de- 
pend on the system variables. If the mean field interactions 
are dominant, the memory function will depend on the system 



8 



wavefunction, which depends on N, This gives the possibility 
of novel linewidth and frequency dependence on the pumping. 



Level shifts and damping rates 

It is worth establishing the connection between our results 
and the more familiar quantum optical theory. 

In terms of the output eigenfunctions, the memory function 
f(t) can be written as 



(61) 



(62) 



f(t) = J dp |(/tK}l e" 1 ^. 
We now evaluate the time integral over f(t), and use 



z + le \ z , 



for e — > 0, where V denotes the principal value part of the 
integral. We find 



du /( u ) e ^o+A SJW )« 

dp \ (k\u p )\ 2 S(u p — u>„ - A sm ) 

K«K)I 2 



-iV / dp 



Up — U! a — A SM 



The damping and frequency shift take the form 

K*K>l a 



A SM = jV / dp 



w + A SM — uj p 



(63) 



(64) 



7sm = IK J dp\(n\u p )\ 5{uj + A SM - uj p ). (65) 

The level shift now takes the familiar form of a principal value 
integral modified by the fact that it is now an implicit equa- 
tion that must be solved simultaneously with (|65|). When 
A SJ/ -C ui , the shift may be found by neglecting A SM in 
the denominator, and can be more easily evaluated using ( |58| ) 
in the same approximation. 



OUTPUT COUPLING 

The potentials in the output field may be caused by mean 
field interactions, gravity, or by a magnetic field. In this sec- 
tion we discuss the memory functions for these process, which 
describe the influence of the output beam on the trapped mode 
statistics. While our model is only strictly valid in the nonin- 
teracting limit where the trapped mode becomes the ground 
state of the harmonic potential, we will also find the memory 
function for the case of a more highly occupied Bose-Einstein 
condensate, because deriving a simple description of this out- 
put coupling process for the Thomas-Fermi regine is of gen- 
eral interest and some practical use. 



Gravity 

If the Greens function for the output potential is known, and 
the position space coupling function (x\k, t) is simple enough 
to compute the necessary integrals, the memory function may 
be found exactly. The memory functions for gravity and the 
magnetically anti-trapped output state have been found using 
the Gaussian form for the coupling [|l7|] 



k(x) 



1 



1/4 



-x 2 /2a 2 



(66) 



When the potential is V(x) = —mgx the result is 



/(*) 



(67) 



where the kinetic and potential energies assiciated with the 
initial wavepacket are hujk = Ti 2 /2mer 2 , and TiLo g = mga 
respectively. The general effect of the output coupling is 
to change the system frequency and the effective damping. 
To show this for the case of gravitational damping we use 
7 = 2 x 10 5 s -2 , and find the Markov and self consistent 
solutions as the pumping is increased. We use n s = 47, and 
increase the pumping from 2 x 10 5 s" 1 to 2 x 10 6 s _1 . We 
use an atomic mass of 5 x 10 _26 kg. The bare system fre- 
quency is taken as uj = 2ir x 123 Hz. We choose gravity 
as g = 9.8 sin(0.1)ms -2 . We will also take the width of the 
coupling to be a - = 1.6 x 10~ 6 m, which corresponds to the 
momentum width used in [O]. Figure || shows the self consis- 
tent steady state occupation and linewdith calculated for these 
parameters. The Markov approximation neglects the effective 
weakening of the damping that is caused by a large number of 
virtual transitions during the output coupling. The more ac- 
curate self consistent approach predicts an increase in TV over 
the Markov results. The difference in linewidths is a conse- 
quence of this effective widening of the system boundary, and 
arises principally from the difference in 7Y. 



Comparison with previous work 

In Figure || we have used parameters corresponding to Ta- 
ble I of pTTj], which examined linewidth as a function of pump- 
ing. The previous treatment used the pump term r instead of 
rN/(N + n s ) in (p4|). When stationary solutions are found, 
instead of the self consistent Markov result (|60|), the linewidth 
becomes 



2N 2(N + n s ) + l 



r(l + 2n s ) 
AN 2 



+ 0{N- 3 ), (68) 



as can easily be verified using the data in Table I. of [|11]]. This 
expression can also be written as 

T SM {l + 2n s ) + 0(N- 3 ), (69) 



2N 2(N + n s ) + l 
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FIG. 2: Numerical results for the Markov approximation (o) and the 
self consistent approach (x). We use 7 = 2 x 10 5 s~ 2 , n s — 47, 
m = 5 x 10 _26 kg. The bare system frequency is taken as uj = 2tt x 
123 Hz. Gravity is set to g = 9.8 sin(0.1)ms" 2 . The width of the 
coupling is a = 1.6x 10 _6 m. a) The steady state occupation number 
as a function of r. b) The self consistent linewidth as a function of r. 



which clearly generates a spurious broadening of the 
linewidth. This masked the physical behaviour found in the 
present work, where the linewidth is reduced below the Born- 
Markov prediction by the non-Markovian evolution. 



Mean field interactions and antitrapping 

When the output atoms are coupled into an antitrapped state 
the memory function should be calculated for an inverted 
parabolic potential. If the atoms also experience significant 
mean field repulsion from a highly occupied trap mode, the 
effective output potential will be described by a rescaled in- 
verted parabola corresponding to the Thomas-Fermi solution 
for the condensate wavefunction. Using the same Gaussian 
form of the coupling and taking the height of the potential at 
the center of the trap as V — Hu, the memory function is [|l7|] 



/(*) 



cosh Lu a t + i 



sinricj £ 

LU LUk 



-1/2 



(70) 



When u 1 « t, this becomes 



/(*) 



-iujt—uj t/2 



Up 



-1/2 



(71) 



from which it is clear that the timescale of decay caused by 
the kinetic dispersion associated with the curvature of the po- 
tential is 2/uj . As noted in |17|, this exponential envelope 



is evident from the Greens function for the inverted parabolic 
potential, and occurs regardless of the shape of the interaction 
region. This is also expected to be the case for the Thomas- 
Fermi mean field potential because while the exact motion 
will not be correctly reconstructed by this potential outside 
the interaction region, it is accurately reproduced where the 
overlap is calculated. 

Finding the output evolution for a particular potential is 
equivalent to finding the Greens function for the output field 
Schrodinger equation. This provides rather more information 
than is needed, since we only require the evolution within the 
region of overlap with the trap state. To find this, we will write 
the output field effective Schrodinger equation in terms of am- 
plitude and phase variables. We may then use the Raman- 
Nath approximation to evaluate the short time evolution when 
the trap is highly occupied since the phase evolution gener- 
ated by the mean field will be dominant. We will also find 
the Raman-Nath time, and the region of validity may then be 
found by comparing this with the memory time arising from 
the approximate memory function. We will first write the out- 
put wave function as ^(x, t) = A(x, i)e lS ^ x '*\ whereby the 
Schrodinger equation reduces to 



dA 
~dt 

dS_ 
~dt 



— (2VS* • VA + AV 2 S) 
2m v ' 



h 

2m 



A 



(VS)' 



7(x) 



(72) 
(73) 



To obtain the effective one dimensional equation of motion, 
the radial wavefunction is assumed to be in the ground state of 
the harmonic trapping potential, with frequency u> r . In terms 
of the single particle condensate wavefunction ip c , the one di- 
mensional Schrodinger equation for the output atoms is 



h 2 d 2 



where the effective potential is 

V(z) = u 12 N\^ c \ 2 = e M 



(74) 

(75) 
(76) 



where this is positive, and zero otherwise. The rescaled in- 
teraction strength is U12 = u\2/2ttR 2 , and the width of 
the ground state is determined by the harmonic oscillator ra- 
dius R = h/mw r . The axial Thomas-Fermi radius is L = 



10 



y/2[i TF /rmjj% and the ratio of trapped-trapped to trapped- 
untrapped s-wave scattering lengths is e, which we take as 
unity in what follows. The effective interaction parameter is 
rescaled by the transverse cross section, and has dimension 
length x energy as required for a one dimensional model. The 
chemical potential becomes 



3Nu u 



mujf, 



2/3 



(77) 



Starting from the density-phase equations for the output 
evolution and reducing to one dimension, we neglect the cur- 
vature term V 2 A/A, and take the initial phase as constant in 
position. The phase equation then depends only on the poten- 
tial, and may be integrated to derive the memory function 



4 V ^cX 



where 



h(x) 



(78) 



(79) 



1 



and ijj c = V c /h. The long time behavior for h(x] 
has been used to factor out the asymptotic form of f(t). 

The Raman-Nath time for the validity of this approxima- 
tion is the timescale of variation of the amplitude A. Using 
only the potential for the phase evolution, and again using the 
smoothly varying envelope approximation so that we drop the 
V5 • VA term, we find the short time solution 



A(z, t) = A(z, 0) exp -cj 



2mL 2 



(80) 



Defining the kinetic energy associated with the length scale of 



the Thomas-Fermi profile as Took 
Nath time is 

1 



h /2mL , the Raman- 



2 



(81) 



This is just the decay caused by the kinetic energy arising from 
the curvature of the potential, as seen in our previous discus- 
sion for the anti-trapped Gaussian wavepacket. Note that this 
does not depend on N because the Thomas-Fermi solution al- 
ways assumes the form of an inverted image of the trap with 
amplitude given by the chemical potential. We can find the 
condition for the validity of the approximate memory func- 
tion for our purposes by comparing this timescale with the 
memory time. The memory time may be defined by the ratio 



R = 



(82) 



When T m is chosen so that the ratio R is small, the contri- 
bution of any time integral equation involving f(t) will be 
significant for times less than T m . We want to find T m when 
R ! <C 1 for the Thomas-Fermi memory function (|78[). We have 



R 



•erf(V-«w c T m ) 



yJ-iuj c T n 



(83) 



We now assume 1 <C u c T m , and use the asymptotic form of 
the memory function [ 20 1 to find 



Auj r R 2 



< T 



Using R = 10 2 , the condition T m <C t rn becomes 



LU < 10 6 UJ C 



(84) 



(85) 



When this condition holds the approximate memory function 
will determine the effect of the output coupling on the system 
evolution. 



Kinetic evolution 

We are primarily interested in the regime where the phase 
evolution determines the memory function. However, if we 
include the long-time exponential decay determined by the 
kinetic evolution, the analysis is greatly simplified. This 
is because neglecting the kinetic evolution entirely leads to 
an unphysical long-time behavior of the memory function, 
which renders the equations of motion insoluble. This fea- 
ture is related to the existence of a bound state in the dressed 
states of the trap plus output potentials, which has been dis- 
cussed in detail elsewhere [[TT|]. Since there is a separation 
of timescales, we may include the kinetic evolution by sim- 
ply multiplying the approximate memory function by e ~ UJ ° t / 2 . 
This has no effect on the overlap /(<) over the region in which 
the integrand is significant, but restores the description of the 
long time behavior to a more realistic form. The memory 
function that we will use takes the final form 



U At) 




h{u) c t)e 



-iuj c t—uj t/2 



(86) 



Mean field damping in the Markov approximation 

We will now find the linewidth and the frequency shift of 
the atom laser model using f TF (t)in the Markov approxima- 
tion, and show that mean field damping leads to an interesting 
scaling behavior for the trap. We can compute an analytic so- 
lution by using the integral 
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dt Mt)e -« = » ( <(' + "./2) ^ C^^-i . (87) 



This is easily found from the integral representation 

UAt) = 11 dx (1 - x^e-^Ki-, 2 ) e -u t/2 > (8g) 
1 Jo 

which compactly expresses the approximations we are us- 
ing for the output evolution. 

We now use the system frequency u = uj c for the con- 
densate, and find that when uj <C w c the asymptotic form of 
arctan [EG] leads to 



37T7 



The occupation becomes 



N = V c - - n s 

37T7 



(89) 



(90) 



so that with the N dependence of uj c = /j, TF /h, well above 
threshold the trap number becomes 



where 



N = C\ - 

,7 



3/2 



3/2 



, _2_ /woV" / 3un Jmuil 



3n y h 



1/2 



(91) 



(92) 



This is a scaling behaviour quite different to the optical case 
which is varies as N ~ r/7. The trap spectrum linewidth is 



CONCLUSIONS 

We have carried out a non-Markovian steady state analysis 
of a fully quantum mechanical atom laser model. We have 
demonstrated that a self-consistent Markov approximation is 
valid provided the laser is operating in a linewidth narrowing 
regime and the reservoir correlation time is sufficiently short. 
We have shown that the difference between the Born-Markov 
approximation and the more exact treatment is that the system 
frequency, occupation number and effective damping may be 
shifted by the coupling process. This leads to corrections in 
the effective damping rate and the steady state trap occupation 
number. We have found a simple analytical form for the mem- 
ory function for a Bose-Einstein condensate output coupler in 
the Thomas-Fermi regime, which gives rise to nonlinear scal- 
ing of the steady state occupation with pumping rate. This has 
implications for the design of linewidth narrowing quantum 
feedback schemes, and for the stability of atom lasers. 
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(93) 



The possibility of the the damping depending on N leads to 
some interesting new scaling behavior. The expression ([JlJ) 
demonstrates that the Thomas-Fermi profile is not a partic- 
ularly good repelling potential, requiring that if the flux is 
increased, the occupation must adjust by a correspondingly 
larger amount in order to repel the extra atoms rapidly enough 
to recover equilibrium. The linewidth scales as j 3 /r 2 , rather 
than the usual optical scaling of r y 2 /r which has consequences 
for the response of atom laser statistics to fluctuations. 

The validity of the Markov approximation is determined by 
the condition A <C u> c < which leads to the requirement 



2/3 



7T7 



< N. 



(94) 



This forces a restriction on the rate of output coupling relative 
to the interaction strength and the occupation. 
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